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Abstract 



f^ ' A realistic medium- or long-baseline neutrino experiment may suffer from terrestrial matter effects 

Cn . which are likely to contaminate the genuine T-violating asymmetry between v^ — >■ Va and v a — ^ v^ 

^ \ oscillations. With the help of the commutators of lepton mass matrices in matter, we show that this 

kind of contamination is negligible for a variety of experiments provided the neutrino beam energy E 



and the baseline length L satisfy the condition 10 ^(L/km)^(GeV/i?) ^ 1. 
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In the framework of three neutrino flavors, the strength of leptonic CP and T violation in neutrino 
oscillations is measured by a universal and rephasing-invariant parameter J [1] defined through 

Im {yo,i^fi^V*jV^^ = ^ X] ^"/97 X] ^yfc ' (^) 

7 k 

in which V is the 3x3 unitary lepton flavor mixing matrix [2], the Greek subscripts run over {e,fj,,T) 
and the Latin subscripts run over (1,2,3). In terms of the standard parametrization of V [3], we have 
J' = sin 6i2 cos 6i2 sin ^23 cos ^23 sin 6^^ cos^ 9^^ sin 6. The fact that the smallest flavor mixing angle 61^ 
is not that small [4] makes it very hopeful to see appreciable CP- and T-violating effects in the lepton 
sector provided the unknown Dirac phase 5 is not strongly suppressed |5]. The Jarlskog parameter J' is 
directly associated with the T-violating asymmetry between the probabilities of v^ — )• Vg and Up — )• z/„ 
oscillations [6]: 

^T = Pi'^a ^^p)- P{^p ^ ^J = 16:^ Sin ^ sin ^ sin ^ , (2) 

in which (a, /3) = (e,/x), (/U,t) or (t, e), Aj- = m? — m^ with m^ or tti- being the neutrino masses (for 
i,j = 1,2,3), E denotes the neutrino beam energy, and L is the distance between the neutrino source 
and the detector. In a realistic medium- or long-baseline neutrino oscillation experiment, however, the 
terrestrial matter effects may contaminate the measurement and thus affect the determination of J (or 
equivalently, 5). It is generally expected that the T-violating asymmetry between v^ — )• Ua and Vo — ;■ v^ 
oscillations should be less sensitive to matter effects than the CP-violating asymmetry between u^ — )• i/^ 
and V^ — )• Vn oscillations [7J, because matter effects on the CP-conserving parts of Piy^ ~^ ^s) ^.nd 
P{i^a —7- i^a) 3-^6 exactly the same. Therefore, we have the T-violating asymmetry in matter [8j: 

^T = PK ^ ^0) - P{^0 ^ ^a) = 16 J^ Sin %^ Sin ^ sin ^ , (3) 

where A • • = mf — tv?- with m^ or m • being the effective neutrino masses in matter, and a "tilde" always 
indicates the quantities contaminated by matter. It is known that i7A^2^23^3i ~ •^^12^23^31 holds for 
a constant matter profile [9], but under what condition can Ar^ = A^^ hold to a good degree of accuracy? 
The present paper aims to answer this question with the help of the commutators of lepton mass 
matrices developed in Ref. [10]. Of course, someone may trivially argue that ^rp must be approximately 
equal to Ar^ if the terrestrial matter effects are small. But then the question becomes how small is small. 
We shall specify the conditions under which A^^ = Ar^ hold either exactly or to a good degree of accuracy. 
Our work is different from Ref. [10] at least in two aspects: (a) we follow a more concise procedure, in 
which some complication and ambiguity can be avoided, to establish the relationship between A^j- and 
^rp; (b) we figure out the explicit condition to assure Ar^ ~ ^rp, which can be used to judge a realistic 
medium- or long-baseline neutrino oscillation experiment. 

Let us work in the flavor basis where the charged lepton mass matrix M^ is diagonal (i.e., M^ = 
Df^ = Diag{m,g,r?T, , m,^}). Then the neutrino mass matrix M^^ can be written as M^^ = VDj^V' with 
Dj^ = D[ag{mi, 7712, m^} . ^^ matter, we denote the effective neutrino mass matrix as M^, and the effective 
lepton flavor mixing matrix as V. Then we have M^ = VD^V^ with D^ = Diag{m]^, 7712,^,3}. In a way 
similar to Ref. [10], we define the leptonic commutators in vacuum and in matter: 

X^= M^M/,exp(^2ieM^iV4) = Z^^^ exp (^21^^1)2^1^ 

X^ = [M,M;,exp (2ieM,Mt)] = [z)|,exp (2\iVDlV^)] , (4) 



where ^ = L/{4:E) with E being the neutrino beam energy and L being the basehne length. When a 
neutrino beam travels through a medium, it sees a nontrivial refractive index because of the coherent 
forward scattering from the constituents of the medium via the charged current interactions [11]. The 
effective Hamiltonian responsible for the propagation of neutrinos in normal matter can be expressed as 



7i = — VD^V^ = — (vd'^V^ + T 



(5) 



where T^ = Diag{j4, 0, 0} with A = 2\/2GpN^E. Here Gp is the Fermi constant and N^ denotes the 
background density of electrons. Eq. (5) implies a simple relationship between M^^Ml and M^Ml; that is, 
VD^V^ = VD'^V^ + ^m- This relationship will allow us to establish a relation between the determinants 
of Xc and Xc under some conditions, as we shall see later on. 

The exponential of a square matrix Z is given by the power series exp(Z) = /+Z+Z^/2! + Z^/3! + - • • , 
where I denotes the identity matrix. Because M^,Ml has three different eigenvalues, we obtain 



exp hiiVDlV^) = Y, ^ „r ^ = V 



n=0 



The commutator Xc turns out to be 



E 

.n=0 



(2iei?^) 



2\" 



n\ 



y^ = yexp(2i^D^^y^ 



X, 



my 



Dj, V exp {2iiDl) V'\=Y.'-^ [^1 VDI^'V' 



n=0 



It is easy to figure out the explicit expression of [Z?!, VDl"'V^ , which must be traceless: 



DJ,VDl''V^ 



A^^Ym^^V^^V^i A,,^m^y,,K* 



i=l 



A Vm2"y*y. 



1=1 
3 



A^,j;m2"y^,X; 



3 



i=l 



in which A^g = m^ — mi (for a, /3 = e, /i, r). As a result, 



X, 



a/3 = '"a - '"/3 ^^" 

3 3 

A,^, Yl «^P (2i?"^') Ve^V;^ A,, ^ ^^P (2^^"^^) K.K* 

i=l i=l 

3 3 

A^, Y exp (2iem2) K*y^, A^, ^ «^P (^i^"^') ^^^^^r^ 

3 3 

A,, ^ exp (2iem2) K*V;, A,^ ^ «^P (^i?"^') ^^^K 

i=l i=l 

A straightforward calculation leads us to the determinant of X^ as follows: 

3 3 3 

det X^ = 2iA,^A^,A,, E E E ^"^ (V;,^^,KfcKfc^;.K*,) exp [21^ (m? + m] + mf,)] 

i=i j,=i k=i 

/ 3 \ 



j=l 



-4 JA.^A^^A,, [sin (2^Ai2) + sin (2CA23) + sin (2eA3i)] exp 21^ ^ 



m,- 



(6) 



(7) 



(8) 



(9) 



(10) 



j=i 



where J has been defined in Eq. (1). The calculation of detXt is exactly analogous, and the final results 



of det Xc and det Xc are 



detXg = 16 JA,^A^,A,, sin -i|^ sin 



Ai2i^,^A23L^.^A3^^^p/ .^^ 



4S 



4S 



4S 



m. 



3 



det X^ = 16:f A,^A^, A,, sin ^ sin ^ sin ^ exp ( 21^ ^ 

\ i=\ 

Comparing Eq. (11) with Eqs. (2) and (3), we arrive at 



m? 
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detXg = ^TAe^A^,A,,exp [ 2iC^mn 
detX^ = ^TAe/.A^,A,,exp (2iC^m2 j 



(12) 



The next step is to establish a relationship between det Xt and det Xp^ in order to establish a relationship 
between A^ and ^rp. Let us now come back to Eq. (5). 

Given two square matrices Y and Z, the exponential oiY ^ Z can be expressed as the Zassenhaus 
expansion [12]: exp (F + Z) = exp (F) exp (Z) exp (- [F, Z\ /2) exp ([Z, [y, Z]] /3 + [F, [y, Z]] /6) • • • . We 
see that exp (Y -\- Z) = exp (Y) exp (Z) exactly holds if Y and Z are commutable (i.e., YZ = ZY). In 
our case, what we are concerned about is the commutator 



im' 



VD?yKT 



V ' ! -^rti 



2fA 



5]a,iK.,f;, E^^i^«K* 



j=2 



i=2 
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i=2 



(13) 



If the condition 



7m = i A 



EA-iF-F* 



i=2 



+ 



E^^i^^K* 



j=2 



<1 



(14) 



is satisfied, one may argue that the commutator in Eq. (13) is far smaller than the identity matrix and 
thus 

exp f 2i^yL>2yt'\ ^ g^p UiiVDlv'^^ exp (2i^Tj„) 



(15) 

holds as a good approximation. Note that exp (2i^Z)3) = Diag {exp (2i^m^) ,exp (2i^rn-2) ,exp (2i^?TT-f)} 
and exp {2\£_T^ = Diag {exp (2i^A) , 1, 1} hold. Therefore, 

■ 3 3 3 

^exp [2i^ (m? + A)] \Vj J^exp {li^mj) V,,V;, J^exp (21^^^) V,,V:, 



exp h.i^VDlv'^^ ~ 



3 



3 



3 



j;exp[2ieK + ^)]y,*y^, J]exp(2iem2)|i;j2 ^exp (21?^^) y^,y; 

3 3 3 

J]exp[2ieK + ^)]K*y,, J^exp(2iem2)v;*y^, J^exp (21^^^) |y„|2 



.i=l 



j=l 



i=l 



(16) 



can be obtained provided 7^^ <C 1. The commutator Xc is then given by 



X^^ 



A 



3 

Me 2^ 
i=l 
3 

re 2_^ 



A,^^exp(2iem2)i;.I/ 

exp[2ie(m2 + yl)]l/*V;, 
exp [2ie (m^ + A)] K*y„ A,^ J^ exp (2iem2) v;,V,, 



i=l 



3 



A,,^exp(2iem2)i;.y;. 
3 



4=1 



j=l 



(17) 



in the same approximation. This allows us to get 



det X^ ~ exp (2i^^) det X^ = ^rp exp 

where we have used the sum rule [13] 
3 



2\i 



'mf + A 



^Texp(2ie^ 



m? 



(18) 



j=i 



Y,^l = tr (VDlV^ = tr [vdIv^ + tJ =Y^m^ + A, 



(19) 



i=l 



i=l 



which can directly be observed from Eq. (5). Comparing Eq. (18) with the second equality in Eq. (12), 
we immediately obtain ^rp ~ Ar^ under the condition j^ <^ 1. Some discussions are in order. 

• Ar^ = Arp exactly holds if M^^Ml = VD'^V' and T^^ are exactly commutable. The latter possibly 
implies: (a) there is no matter effect, ^ = 0; or (b) the neutrino masses are exactly degenerate, 
nil = m2 = m^; or (c) there is no lepton flavor mixing, V = I. In either case (b) or case (c), of 
course, there would be no neutrino oscillation phenomenon. 

• The order of j^ can be estimated as follows. A global fit of current neutrino oscillation data 
yields IA31I ~ 3OA21, K,V;,\ ~ 3|K3^;3l, 1^2^*21 ~ 31^3^*31, and \V^,\ 



~ \V^^\. Therefore, 



j^^2eA\A,,\\v,,v;,\r.w-'(:^ 



rm 



(20) 



where jAgi] ~ 2.4 x lO"'^ eV% \V^^\ ~ 0.15, |V;,3| ~ 0.7 iteFit, and A ~ 2.28 x IQ-'^eY^E/GeV [E] 
have typically been input. Given the T2K experiment, one has L ~ 295 km and E ~ 0.6 GeV, and 



thus 7j^ 



3. As for the NOz^A experiment, L ~ 810 km and E c^2 GeV, leading us to 7^ 



So 7j^ ^ 1 is satisfied in both of these two experiments. 

If the condition 7^^ ^ 1 is not satisfied, then the commutator in Eq. (13) and the higher order expansion 
terms will contribute an appreciable correction to Eq. (15). In other words, A^^ should contain appreciable 
matter effects in this case, making it difficult to extract clean J ox 5 ui this kind of experiments. 

The above calculation is subject to a neutrino beam traveling in matter, but it can directly be 
extended to the case of an antineutrino beam traveling in matter. Thanks to CPT invariance, we have 
P{pa ~^ ^ r) — ^(^« ~^ ^a) a^*^ ^^ 13 ~^ ^a) = P{''^a ~^ ^b) i^ vacuum. The T-violating asymmetry 
between V^ — )• Va and Va — )• V^ oscillations is therefore given by ^^ = —Ar^, as one may easily read off 
from Eq. (2). In matter, the corresponding T-violating asymmetry A'r^ can be approximately equal to 
A!rY provided the condition 7^^ ^ 1 is also satisfied lj. 



""^As for an antineutrino beam in vacuum or in matter, one may simply make the replacements V -^ V , J" — !■ —J" and 
A — )■ —A in the relevant calculations. 



Note that the CP-violating asymmetry A^p = P{i^a ~^ ^b) ~ P^a ~^ ^b) i^ exactly equal to 
the T-violating asyinnietry Ar^ in vacuum, as guaranteed by CPT invariance. In matter, however, 
Aqp = P{i^a ~^ ^r) ~ ^{^a ~^ ^r) is in general different from ^rp because the CP-conserving parts 
of P{ya ~^ ^r) ^iid -P(^'„ —7- Va) are not identical with each other due to the opposite-sign matter 
corrections q This kind of fake CPT violation is actually a pure matter effect [16] , which may more or 
less contaminate the extraction of the genuine CP-violating effect from ^^p [IZ] • 

Of course, it is always possible to make an analytical expansion of P(f„ — )• Vo) in terms of small 
parameters (e.g., the ratio A2i/A3^ and sin^j^g) [HlIISj, so as to see the sensitivities of its CP-conserving 
(or T-conserving) and CP-violating (or T-violating) parts to terrestrial matter effects. In the present work 
we have provided an alternative way towards understanding matter effects on the T-violating asymmetry 
in neutrino oscillations. This way is certainly more compact, although its usefulness is limited to the 
cases in which 7^^ <C 1 must be satisfied. Fortunately, it seems that most currently-proposed medium- 
and long-baseline neutrino oscillation experiments are consistent with this condition. 

In summary, the discovery of leptonic CP or T violation will be one of the major targets in tomorrow's 
experimental neutrino physics. To achieve this important goal, a number of long- or medium-baseline 
neutrino oscillation experiments have been proposed or are under construction. We have shown that 
most of such experiments satisfy the condition 10^ (L/km)^(GeV/£') <^ 1, implying that the T-violating 
asymmetry between v^ — )• Va and Va — )• v^ oscillations (or between V^ — )• Va and Vo — )• V^ oscillations) 
is essentially free from terrestrial matter effects. In our proof we have adopted the language of leptonic 
commutators in matter, which can provide us with a concise and novel way of understanding matter 
effects. We hope that this useful theoretical language could find some more phenomenological applications 
in neutrino physics. 

I am grateful to Y.F. Li and Y.L. Zhou for stimulating and helpful discussions, and to T. Ohlsson for 
interesting communications. This work is supported in part by the National Natural Science Foundation 
of China under Grant No. 11135009. 



^Namely, the effective Hamiltonian responsible for the propagation of antineutrinos in normal matter is given by T-L'^ 
[VDlV'^ - T^) I {2E), in contrast to U^ = {VDlV'^ + T^) / {2E) for neutrinos. 
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